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HD 80606 b, a planet on an extremely elongated orbit^ 
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Abstract. We report the detection of a planetary companion orbiting the solar-type star HD 80606, the brighter 
component of a wide binary with a projected separation of about 2000 AU. Using high-signal spectroscopic 
observations of the two components of the visual binary, we show that they are nearly identical. The planet has 
an orbital period of 111.8 days and a minimum mass of 3.9Mj up . With e= 0.927, this planet has the highest 
orbital eccentricity among the extrasolar planets detected so far. We finally list several processes this extreme 
eccentricity could result from. 
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1. Introduction 

We report in this paper on our radial-velocity measure- 
ments of HD 80606, the primary star of the visual binary 
system HD 80606-HD 80607. These observations reveal the 
presence of a 3.9 Jovian-mass planet (minimum mass) in 
, a very eccentric orbit around this solar-type star. 

The variable velocity of HD 80606 was fir st noticed 
by the G-Dwarf Planet Search ( Latham 200C ) , a recon- 
naissance of nearly 1000 nearby G d warfs that uses the 
HIRES high-resolution spectrograph ( Vogt et al. 1994 ) 



mounted on the 10-m Keck 1 telescope at the W.M. Keck 
Observatory (Hawaii, USA) to identify extrasolar planet 
candidates. The star was then followed up by the ELODIE 



Planet Search Survey team (Mayor & Queloz 1996; Udry 
et al. 2000) using the ELODIE fiber-fed echelle spectro 
graph (Baranne et al. 1996) mounted on the Cassegrain fo 
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cus of the 1.93-m telescope at the Observatoire de Haute- 
Provence (CNRS, France). 

The ELODIE velocities are obtained by cross- 
correlating the observed spectra with a numerical tem- 
plate. The instrumental drifts are monitored and corrected 
using the "simultaneous Thorium- Argon technique" with 
dual fibers ( Baranne et al. 1996 ). The achieved precision 
with this instrument is of the order of 10 m s~ 1 . The HIRES 
instrumental profile and drifts are monitored using an 
Iodine gas absorption cell (Marcy & Butler 1992). The 
radial velocities are derived from the spectra using the 
TODCOR code ( pucker fc Mazeh 1994] ), a two-dimensional 
correlation algorithm. 

The observations of HD 80606 started in April 1999 
with HIRES. With a velocity difference of 267ms~ 1 in 
less than one month between the first two measurements, 
the variability of this source was quickly detected. In 
July 1999, we started an ELODIE radial-velocity follow 
up of 6 non-active slow-rotating radial-velocity variable 
stars detected with HIRES, including HD 80606. The first 
ELODIE measurement for this star was obtained dur- 
ing our November 1999 run. The discovery of the plan- 
etary companion orbiting HD 80606 has been recently an- 
nounced together with 10 other new extrasolar planet 
candidates (April 4th 2001 ESO PR0). Among these is 
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the planetary companion to HD 178911 B (Zucker et al. 
in prep.), another one of the candidates identified with 
HIRES. 

The stellar characteristics of the two components of the 
HD 80606-HD 80607 visual binary are presented in Sect.||. 
The radial-velocity data and the orbital solution are pre- 
sented in Sect.H. The very high orbital eccentricity is dis- 
cussed in Sect. ET The 61 radial-velocity measurements pre- 
sented in SectTH as well as the Iron line list we used in 
Sect. ^| will be made available in electronic form at the 
CDS in Strasbourg. 

2. Stellar properties 

HD 80606 (HIP 45982) and HD 80607 (HIP 45983) are the 
two components of a visual binary system. They have 
common proper motions and the fitted systemic ve- 
locity for HD 80606 (7 = 3.767± O-OlOkms" 1 ) is almost 
equal to the mean radial velocity measured for HD 80607 
((RV) =3.438±0.025kms" 1 ). The difference between the 



Table 1. Observed and inferred stellar parameters for 
HD 80606 and HD 80607 



HD 80606 



HD 80607 



Sp. Type 








G5 






G5 


mv 






9.06 


± 


0.04 


9.17 


± 0.04 


B — V 






0.765 


± 


0.025 


0.828 


± 0.029 


7T 


(mas] 




17.13 


± 


5.77 


9.51 


± 8.76 


Distance 


(pc) 




58.4 


± 


29.6 
14.7 


105 


1 1228 
± 51 


cos (5) 


(mas 


yr" 1 ) 


46.98 


± 


6.32 


42.90 


± 9.23 




(mas 


yr" 1 ) 


6.92 


± 


3.99 


8.26 


± 5.88 




(°K) 




5645 


± 


45 


5555 


± 45 


log5 


(cgs) 




4.50 


± 


0.20 


4.52 


± 0.15 




(km s 




0.81 


± 


0.12 


0.91 


± 0.12 


[Fe/H] 






0.43 


± 


0.06 


0.38 


± 0.06 


vsini 


(km s 




0.9 


± 


0.6 


1.4 


± 0.4 


W XM 


(mA) 




<2.5 


<3.0 


logn(Li) 






<0.78 


<0.77 



two values can be explained by the binary orbital mo- 
tion. The main stellar characteristics of HD 80606 and 
HD 80607 are listed in Table||. The spectral types, appar- 
ent magnitudes, colour indexes, parallaxes and proper mo- 
tions are from the HIPPARCOS Catalogue flESA 1997| ). 
The projected stellar rotational velocity, i>sini, was mea- 
sured using the mean ELODIE cross-correlation dip width 
and the calibration by Queloz et al. (19~98| ). The rms of 



widths for both stars. The abundance upper limits were 
then derived using the curves of growth by Soderblom 
(1995 ). The Lithium abundances are scaled with 



et al. 



the HIPPARCOS photometric data is large for both stars 
(<th p — 40 mmag) but this measured scatter is classified 
as 'duplicity-induced-variability' in this catalogue. The 
angular separation between the two visual components is 
about 30 " . This value is not much larger than the satellite 
detector size so contamination from one component onto 
the other is probably responsible for the observed scatter. 
The contamination is also probably responsible for the dif- 
ference in parallaxes (a factor of two) and for the abnor- 
mally large uncertainties on this parameter (av ~ 1 mas is 
expected with HIPPARCOS for a 9th magnitude star). 

We derived the atmospheric parameters (LTE analy- 
sis) using HIRES hi gh signal-to-noise sp ectra with the 
same method as in (Santos et al. (2000a|) . We used the 
same line list and oscillator strengths as these authors, 
except for some lines that could not be used because they 
were out of the HIRES spectral coverage or fell just be- 
tween two non-overlapping orders of the cchelle spectra. 
Our line list finally consisted of 18 Fei lines and only 3 
Fen lines. We estimated the uncertainties on the derived 



atmospheric parameters in the same way as in Gonzalez 
fc Vanturc (1998] ). The two stars have almost the same 



Iron abundance and are very metal-rich dwarfs (respec- 
tively 2.7 and 2.4 times the solar Iron abundance). An 
independent study (Buchhave et al. in prep.) using the 
same HIRES spectra but a different line list gives consis- 
tent results. For the Lithium abundance measurement, we 
summed all our ELODIE spectra in the A 6707.8 A Lil line 
region. No trace of Lithium was detected giving upper lim- 
its (3-cr confidence level) on the corresponding equivalent 



logn(H) = 12. 

3. Radial-velocity analysis and orbital solution 

On the 24th of April 2001, we had in hand a total of 61 
radial-velocity measurements for analysis: 6 from HIRES 
and 55 from ELODIE. The mean uncertainty on the veloc- 
ities are of the order of 14ms -1 (systematic error + pho- 
ton noise) for both instruments. The HIRES velocities have 
an arbitrary zero point. From contemporaneous observa- 
tions, we applied a preliminary shift to these velocitites to 
bring them to the ELODIE system: ARV = +3.807kms~ 1 . 
To account for possible errors in this zero-order shift, the 
orbital solution presented in Table|| includes the resid- 
ual velocity offset ARV H -e between HIRES and ELODIE 
as an additional free parameter. The obtained ARVh-e 
is consistent with zero. Figure |l|a shows the temporal ve- 
locities for HD 80606. The phase- folded velocities are dis- 
played in Fig[j]c. The fitted orbital eccentricity is ex- 
tremely high — e = 0.927 ±0.012. Assuming a mass of 
1.1 Mq for HD 80606, a typical value for a very metal-rich 
star with a solar effective temperature, the planetary com- 
panion minimum mass is mi = 3.90 ± 0.09 Mj up . The semi- 
major axis is 0.469 AU and the orbital separation ranges 
from 0.034 AU (periastron) to 0.905 AU (apastron). 

The residuals to the fitted orbit cannot be explained 
by our measurement errors. The computed x 2 probability 
for the full set of data is lower than 10~ 3 (x 2 = 102.45, 
v= number of degrees of freedom =N—7 free param- 
eters =54). Our measurement errors are correctly esti- 
mated for both instruments (see e.g. the low residuals 
value obtained for HD 178911 B, Zucker et al. in prep.). 
The very low P(x 2 ) value found for HD 80606 can there- 
fore not result from an underestimation of our measure- 
ment errors. Using our HIRES high-signal spectrum, no 
chromosphcric emission is detected for HD 80606 so the 
|cxpccted stellar jitter is low (a few ms' 1 , see e.g~3anto£ 
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Table 2. Fitted orbital elements to the radial- velocity 
measurements for HD 80606. The velocities obtained with 
the HIRES spectrograph (H) have been set into the 
ELODIE (E) system 



p 


days 


111.81 ± 0.23 




T 


HJD 


2 451973.72 ± 0.29 


> 
K 


e 




0.927 ± 0.012 




7 


kms -1 


3.767 ± 0.010 




w 





291.0 ± 6.7 




Ki 


ms _1 


411 ± 31 




ARVh-b 


ms" 1 


1.5 ± 8.5 


j: 


ai sini 


10~ 3 AU 


1.581 ± 0.037 


o 


AM ^ 


1O" 8 M 


4.26 ± 0.29 


i 

o 


7Ti2 sin i 


Mj up 


3.90 ± 0.09 




N 




55(E) + 6(H) 




0"O-C 


ms -1 


17.7 (E:16.3, H:29.9) 





et al. 2000b; 3aar et al. 1998). Activity related processes 
are therefore probably not responsible for the observed 
residuals. The later could be explained by the presence 
of another planet around HD 80606 on a longer period or- 
bit perturbating the stellar radial-velocity signal induced 
by the inner companion. No clear velocity trend was de- 
tected from the residuals curve (see Fig.jlfb). Future mea- 
surements should help to solve the question. 

4. Discussion 

The fitted orbital eccentricity is the highest found so 
far for an extrasolar planet orbiting a solar-type star. 
The orbital eccentricities for extrasolar planets with pe- 
riod lo nger than 100 days a l most cover the full p ossible 
range flMayor fc Udry 2000| ; |Udry fc Mayor 200l|) : from 
nearly circular (see e.g. the recently announced planet 
around HD28185, P = 385 days, e = 0.06, ESO PR 2 ) to 
nearly unity, as in the case of HD 80606. The distribution 
of the eccentricities of the planetary orbits might be a key- 
stone in understanding the formation processes of planets, 
as was pointed out early in the study of extrasolar planets 
by [Mazeh et al. (1997b[ ). 

Before discussing any mechanism that could have gen- 
erated the eccentricity of HD 80606, it is interesting to 
note that the eccentricity distribution of the planets with 
long orbital periods found so far is strikingly similar to 



that of the binary orbits ( 


Heacox 199E ; 3tepinski & Black 


200C, 


2001; 


Mayor & Udry 2000|; Mazeh & Zucker 2000]). 



In particular, the high eccentricity of HD 80606 is very 
similar 



4.4 

4.2 
4 
3.8 
3.6 

40 




_L 



250 



500 750 
HJD-2 451 000 



1000 




Lo one of Lhc highest eccentricity found so fat 1 for 



a spectroscopic binary — 0.975 ( Duquomioy it al. 1992| ) 



Fig. 1. HD 80606 radial- velocity data. Crosses: Elodie- 
OHP measurements. Open squares: Hires-KECK mea- 
surements, a) Temporal velocities, b) Residuals around 
the solution, c) Phase-folded velocities 



2001) strongly suggests that we are dealing with two dis- 
tinct populations. Nevertheless, we might need to look for 
mechanism(s) that can produce a range of eccentricities 
from zero up to unity for the two populations. 

A mechanism to generate eccentric orbits could be the 
gravitational interaction of a planet (and a binary) with a 
disk (Artymowicz et al. 1991; Artymowicz 1992). However, 
a recent study ( Papaloizou et al. 200l| ) suggests that for 
a standard disk model this can happen only for massive 
companions, at least in the range of brown dwarf masses. 
For companions with planetary masses the disk probably 
acts to damp the eccentricity growth, and therefore can 
not explain the observed high eccentricities. 

Another possible mechanism is the gravitational inter- 
action with another planet (s). This could be via dynamical 
instability (e.g. Wcidcnschilling fc Marzari 199(: ; Rasio & 
Ford 1996| ; [Lin fc Ida 1997| ; |Ford et al. 2001|) or through 



The similarity of the two eccentricity distributions does 
not prove that the planets and the low-mass stellar com- 
panions come from the same population. The large gap 
between the mass distrib ution of the planets and that 
of the stella r companions ( Jorissen et al. 2001 ; Zucker fc 
Mazeh 2001 ) and the differences in the met allicity distri- 
butions for stars with and without planets ( pantos et al 



| 2 | www. eso.org/ outreach / press-rel / pr-2001 / pr-07-01 .html 



s ome resonant intera ction with a disk and another planet 
( Murray et al. 2001 ). The instabilities naturally lead to 
high eccentricities, specially if they involve ejection of an- 
other planet out of the system. The resonant interaction, 
on the other hand, seems to need some fine tuning for 
generating eccentricities as high as the one found here. 

A possible clue to the origin of the particularly high 
eccentricity found here could have been found in the fact 
that HD 80606 resides in a stellar wide binary. At least one 
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other planet, the one orbiting 16CygB, was found with a 



high eccentricity in a wide binary. A few studies (Mazeh 



et al. 1997a; Holman et al. 1997) have suggested that the 
high eccentricity of 16CygB is because of the gravita- 
tional interaction with the distant stellar companion. In 
this model, the 'tidal' interaction of the distant compan- 
ion induces an eccentricity modulation into the planetary 
orbit on a long timescale. The present phase of the cycle is 
close to the highest point of the eccentricity modulation. 

However, this interpretation does not apply here. 
This is so because the modulation ti mescale induced by 
HD 80607 is of the order of 1 Gyr (e.g. |Mazeh & Shaham 



1979()~frhis is very long relative to the relativistic perias 
tron passage modulation, which is of the order of 1 Myr. 
The precession of the longitude of the periastron induced 
by the relativistic effect completely suppresses the third- 
body modulation. To hold on to the third-body interpre- 
tation we have to assume an additional body in the sys- 
tem, in an orbit around HD 80606 with a period of the or- 
der of 100 yrs. The present radial-velocity measurements 
can not rule out such a companion. To be consistent, this 
model has to apply for all the planets with high eccentric- 
ity, above, say, 0.6 — an eccentricity that does not seem to 
be that rare anymore. Note also that this model requires 
a large angle between the plane of motion of the planet 
and that of the perturbating body. 

In short, as stressed by Mayor et al. (2000) and 
Stepinski fe Black (200l| ), we need a consistent model 



that will account for the distribution of eccentricities of 
the planetary orbits and for its similarity to the distribu- 
tion for stellar companions. It seems that further observa- 
tions and theoretical work are needed to reach a consensus 
about such a model. 
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